State reconstruction for composite systems of two spatial qubits 
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Pure entangled states of two spatial qudits have been produced by using the momentum transverse 
correlation of the parametric down-converted photons [Phys. Rev. Lett. 94 100501]. Here we show 
a generalization of this process to enable the creation of mixed states of spatial qudits and by using 
the technique proposed we generate mixed states of spatial qubits. We also report how the process of 
quantum tomography is experimentally implemented to characterize these states. This tomographic 
reconstruction is based on the free evolution of spatial qubits, coincidence detection and a filtering 
process. The reconstruction method can be generalized for the case of two spatial qudits. 
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I. INTRODUCTION 

The concept of quantum state plays a central role in 
the Quantum Theory. It is considered to be the most 
complete description for a physical system. Statistical 
distributions of results of experiments carried out on a 
physical system can be completely predicted from its ini- 
tial state. Therefore, the experimental determination of 
an initially unknown quantum state becomes a very im- 
portant research subject. This has led to the develop- 
ment of techniques to perform the state determination. 
In recent years the problem of state determination has 
received a considerable degree of attention due to impor- 
tant results in Quantum Information Theory. 

Several techniques have been designed and used for 
the state estimation of different physical systems. In the 
field of atomic physics, quantum endoscopy was used to 
determine the state of ions and atoms [H, 0, Q • In quan- 
tum optics, the Wigner function of multi mode fields 
could be measured using homodyne detection 0, d, 
and the technique of quantum tomographic reconstruc- 
tion (QTR) was used for measuring the polarization state 
of parametric down-converted photons [7|. 

In general, these methods are based on a linear inver- 
sion of the measured data. In the case of QTR, the data 
is acquired with a series of measurements performed on 
a large number of identically prepared copies of a quan- 
tum system. The fact that this transformation is linear, 
makes it strongly dependent of any experimental error 
that may occur while recording the data. It can appear as 
a consequence of the experimental noise or misalignment 
and therefore, the reconstructed state is only a reasonable 
approximation of the real quantum state. The density 
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matrices obtained may have properties that are not fully 
compatible with a quantum state. Another alternative 
that has been considered for the state determination is 
the numerical technique called maximum likelihood es- 
timation [H, It is based on a relation between the 
measured data and the quantum state that could have 
generated them. Even though it generates only possible 
density matrices, it has the drawback of enhancing the 
uncertainty on the state estimation. 

In this article we are interested in the determination of 
the state of a composite system. In our experiment, the 
system corresponds to two photons generated by sponta- 
neous parametric down-conversion (SPDC). In this non- 
linear process a photon from a pump laser beam incident 
to a non-linear crystal originates probabilistically two 
photons, signal and idler [1 01 ] . The photons of this pair 
are also called twin photons for being generated simul- 
taneously [ll|. Recently, we have demonstrated that by 
placing D symmetric slits at the path of each twin photon 
it is possible to use the transverse correlations of the pho- 
ton pair to generate maximally entangled states of two 
effective D-dimensional quantum systems [HI, [l3|, fl4j . 
We refer to these D-dimensional quantum systems as 
spatial qudits. In the present work we extend previous 
results to the generation of mixed states of two spatial 
qubits, which also applies to the case of qudits. Following 
this, we investigate the state determination of two spatial 
qubits. We show both theoretically as well as experimen- 
tally, that one can implement the process of QTR to ob- 
tain the density operator of a state composed of two spa- 
tial qubits. The quality of the reconstruction performed 
is also discussed. Even though we had considered only 
the special case of spatial qubits, it is straightforward to 
show that the technique used can be generalized for be- 
ing applied to a system composed of two spatial qudits. 
The main motivation on studying both the generation 
and the reconstruction of mixed states of spatial qudits 
is to consider more realistic experimental situations in 
case of using them in technological fields, such as quan- 
tum communication, where pure states can become into 
mixed ones due to interactions with their environment. 
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II. CONTROLLED GENERATION OF MIXED 
STATES 

It was shown in [T3, EH that the state of parametric 
down-converted photons when each photon is transmit- 
ted through identical multi slits is given by 
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where D is the number of slits in each multi slits array, 
d is the distance between the center of two consecutive 
slits, a is the half width of the slits and Id = {D — l)/2. 
The function Wi m is the spatial distribution of the pump 
beam at the plane of the multi slits (z = za) and at the 
transverse position x — (I + m)d/2, 



W, m = W 



(I + m)d 



-',za 
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The \l) 8 (or \m)i) state is a single-photon state de- 
fined, up to a global phase factor, by the expression 



0, 
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and represents the photon in mode j (j = i,s) trans- 
mitted by the slit The transverse component of wave 
vector of the down-converted photons in the mode j is 
represented by qj . The states in the set { | / )^ } are or- 
thonormal, that is — 8w . We use these states 

to define the logical states of the qudits. In this sense 
(H|) describes a composite system of two qudits. Each 
qudit is represented by a state in a Hilbert space of di- 
mension D, being D the number of available paths for its 
transmission through the multi slits array. 

It can be seen from ([T]) and , that it is possible to 
create different pure states of spatial qudits if one knows 
how to manipulate the pump beam in order to generate 
distinct transverse profiles at the plane of the multi slits 
(W(£; za)) [HI- Let us now assume that, before reaching 
the crystal, the pump beam pass through an unbalanced 
Mach-Zehnder interferometer where the transverse pro- 
file of the laser beam is modified differently in each arm. 
If the difference between the lengths of these arms is set 
larger than the laser coherence length, we will obtain an 
incoherent superposition of the spatial qudits states gen- 
erated by each arm. 

We show in the following section how to use the QTR 
technique to determine the density matrix of these com- 
posite systems. The state whose the density matrix is 
reconstructed experimentally is a mixed state of two spa- 
tial qubits. This state is generated with the experimental 
setup represented in figure flja). A 5 mm /3-barium bo- 
rate crystal is pumped by a 500 mW krypton laser emit- 
ting at A = 413 nm for generating SPDC. Before being 
incident at the crystal, the pump beam cross an unbal- 
anced Mach-Zehnder interferometer. The difference be- 
tween the lengths of each interferometer arm (200 mm) is 



set larger than the laser coherence length (80 mm) . Two 
identical double slits A s and A, are aligned in the direc- 
tion of the signal and idler beams, respectively, at a dis- 
tance of 200 mm from the crystal {za)- The slits' width 
is 2a — 0.09 mm and their separation, d = 0.18 mm. 
The smaller dimension of the double-slits are in the x- 
direction. All measurements are done in the x-axes, at 
the detection plane. At the arm 1 of the interferometer, 
we place a lens that focus the laser beam at the plane of 
these double slits, into a region smaller than d. In arm 2, 
we use a set of lenses that increases the transverse width 
of the laser beam at za- The transverse profiles gener- 
ated are illustrated in figure QTa). The photons transmit- 
ted through the double-slits are detected in coincidence 
between the detectors Di and D s . Two identical sin- 
gle slits of dimension 5.0 x 0.1 mm and two interference 
filters centered at 826 nm with 8 nm full width at half 
maximum (FWHM) bandwidth are placed in front of the 
detectors. 
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FIG. 1: (a) Schematic diagram of the experimental setup used 
for generating and characterizing the mixed states of spatial 
qubits. The pump beam that cross arm 1 is focused in a 
narrow region at za or in a broader spatial region when it 
cross arm 2. A s and Ai are the double-slits at signal and idler 
propagation paths, respectively. D 3 and Di are detectors and 
C is a photon coincidence counter. The configuration used to 
determine the diagonal elements is represent in (b). (c) and 
(d) were used for the second type of measurement and (e) for 
the third type. All measurements are done in the x-axes, at 
the detection plane. 

By using {T]) and ([2|) , we can show that the two-photon 
state, after the double slits, when only arm 1 is open is 
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given by 

l*)i = ^(l + ).|-)i + |-)-l+)0- (4) 

To simplify, we used the state | + )j and | — )j in place of 
the states | \ )j and | )j- Thus, the state (| — )j) 
represents the photon in mode j being transmitted by 
the upper (lower) slit of its double slit. The state of ^ 
is a maximally entangled state of two spatial qubits and 
their correlation is such that when the idler photon passes 
through the upper (lower) slit of its double slit, the signal 
photon will only pass through its lower (upper) slit, and 
vice- versa. 

However, if the laser beam cross only arm 2, the state 
of the twin photons transmitted by these apertures will 
be given by 

l*> 2 = ie < *(|-).|+) i + | + ).|-) i ) 

+5(|->-|->i + l+>-l+>i), (5) 

where <j> = kd 2 /8zA- And now the correlation of this 
spatial qubits is different since we can also have both 
photons of the pair generated, crossing their upper or 
lower slits, simultaneously. 

Therefore, the two-photon state generated in our ex- 
periment when the two arms are liberated, is a mixed 
state of the spatial maximally entangled state shown in 
(j4|) and the state of d5j) . It is described by the density 
operator 

p the = A\ * > X1 <*|+ B\ *> 22 <*|. (6) 

where A and B are the probabilities for generating the 
states of arm 1 and arm 2, respectively. 

III. RECONSTRUCTION 

Now, wc show how QTR can be experimentally imple- 
mented to reconstruct the density operator of the state 
^ generated in our setup without the use of any infor- 
mation about the scheme used for this generation. 

Let us briefly review the process of quantum tomog- 
raphy. The diagonal elements of any density operator 
can be measured directly. Therefore, quantum tomo- 
graphic reconstruction is a protocol to determine the 
non-diagonal elements. It consists in the use of known 
unitary transformations on the system. Each transfor- 
mation generates a new density operator whose diago- 
nal elements are a combination of the coefficients of the 
transformation and of the non-diagonal elements of the 
original density operator. These new diagonal elements 
can be measured. The iteration of this procedure creates 
a set of equations which allows the determination of the 
non-diagonal elements of the initially unknown density 
operator. For a detailed account on this subject we refer 
the reader to [lj| [l6[ ■ 



To characterize the density operator generated in our 
experiment, we first adopt a general form for it 

P = X! PisUm smi \ls,k) (m s ,mi\ . (7) 

l,m— ± 

As we are dealing with two qubits states, the total 
number of measurement basis necessary for the QTR is 
nine [171 ] . They can be generated by using three basis 
for each qubit. Since we have spin 1/2 like systems, 
these basis are the eigenvectors of the Pauli operators 
{cr x ,cr v ,cr z y. In our case, the eigenvectors of a z are the 
slit's states given by ([3]). They form the logical base, 
such that Uj Z ^ = Ij. Measurements in basis {a x ,cr y } 

in mode j require local unitary operations {U^ x \ U^} , 
which allow us for going from a z eigenvectors to a x , a y 
eigenvectors, respectively. The density operator under 
these local transformations can be written as 

with X,fi = x, y. Since the a z eigenvectors are the slit's 
states, the question which remains is: How one can im- 
plement the above discrete local operations in these slit's 
states to perform the QTR of pi For answering this, we 
first consider the state of a photon crossing a given slit I 
(I = ±) along mode j (j = s,i), and propagating through 
the free space to a detection plane located at position z. 
This state can be calculated by the method presented in 
[20l | and is described by 



= y- j d 1j exp(-iaq?)exp(-ilqjd) sine (qja)\ lq 3 ■). 

(9) 

It corresponds to the free evolution of the state \l)j gen- 
erated by the unitary operator Uj (restricted to the one- 
photon subspace), that is \gi)j = Uj\l)j. The operator 
Uj is given by 

Uj = exp(— ik(z — Za)) Jdq exp(~iaq 2 )\lq)jj(lq\, 

(10) 

where a = (z — ZA)/2k . 

Therefore, if we propagate the state p from the plane- 
za to plane-z the result is 

Pz = ^2 Pl*h*n,mi\9l.,9h) (9m s ,9mi\, (H) 
l.m— +,— 

and it becomes clear that p and p z have the same coeffi- 
cients, and thus, that one can reconstruct p by determin- 
ing p Zl i.e, by doing the measurements in the detection 
plane-z. 

It can be deduced from © that photons spread out 
along the measurement plane, so that we have passed 
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from discrete variables, states \l)j, to a continuously dis- 
tributed state | gi )j. Thus, for carrying out the measure- 
ment in discrete basis in plane-z, we need to implement 
an adequate postelection process, which will allow us to 
recover the discrete nature of the logical states. As we 
shall show in the following lines, this can be properly 
done by allocating at the detection plane two new slits 
for each mode j. 

In order to explicitly show how the transformations 



(A) 



, and (A, fj, = x,y) can be implemented to do 



U. 

the QTR for p z , we first write an arbitrary two photon 
pure state in the transverse plane- z as 



*) z oc ^ A m ^ n \g m ) s \g n 



(12) 



T(qj) is the Fourier transform of the double slits trans- 
mission function at mode j 



T^qj) oc ( e "»^° + e^*"- 1 ) sine {qjb) , (15) 

where is the position of the slit k (k = 0, 1) in mode 
j at the plane-z. By replacing the expressions for the 
states | <7; )j in | \& ) z , we determine the biphoton ampli- 



tude F(q s ,qi) in l[T2"|). By inserting F(q s ,qi) and (fl~5|) 
into (| 14(1 . we obtain Ft x , (q s , qi)- After a straightforward 
derivation, we can rewrite the two photon transmitted 
state (fT5|) as 



m,n— ± 



The transmitted state through the double slits placed 
at the detection plane is 



oc 



E 

k,l=0,l 



B k ,i\f{x\,k))s\f{x^))i, (16) 



l*0»= dq s d qi F T ^( qs , qi )\l qs )\l qi ), (13) 



where 



where the transmitted state biphoton amplitude is given 
by 



Ft Xm (q a ,Qi) = Jdq' s J H% F(q' a ,q'i)T\(q' a - q s )T^{q' t - qi ) 



B k ,i = ^ r ™( x ><-k)r n {x^i)A rt 



(17) 



The states of the post-selected photons which crossed this 
(14) additional pair of slits are described by 



J 



= J d 9j exp(-igja; /1>fc )sinc (ij^ + 1j b ) ), 



(18) 



where x^k is the position of the slit k (k — 0, 1) in mode j 
at the plane-z. The transversal position of these slits de- 
termines which effective unitary operation was performed 
at the transmitted photon. We have also defined 



r±(x 



exp 



4a 



2a 



(19) 

Here, by comparing (fl6|) and (fl"2|) . it can be observed 
that the post selection process acts on the | g± )j states 
with the following effective transformation 



U^IU)] « r± (z„,o)|/(a?/i,o)> j 



■^±(^,1)1/(^,1))^ 
(20) 

where | <?± ) j state denotes the post selected state arising 
from I g± )j state. By considering the value of the exper- 
imental parameters: z — za, d, a and &, it can be shown 
that the states |/(x Mj0 ))j an d are orthogonal 

when the condition — a^ol > 4& is satisfied. 



The positions of the new slits for generating the ef- 
fective transformation (Uj) are x x $ = [x Vt o = 
-A/2) and x xA = A (x yA = A/2), with A = = 1376 
mm (Note that condition \x^% — x^^\ — A > 46 is widely 
satisfied). We remark that we refer to the transforma- 
tions done by these slits as effective transformations, 
due to the fact that they act as unitary transformations 
only for the photons in mode j which were transmitted 
through the pair of slits at the detection plane. 

Because of the linearity of quantum mechanics and 
because we are performing only local operations to the 
twin photons, we know that the diagonal elements of 
p era tors are linear combinations of the coefficients 
Pi s iim s mi of pip . The diagonal elements of the trans- 
formed density operators pv*>W a t the plane-z are simply 
determined by using four coincidence numbers measured 
when D s is at the position x\.k for k = 0, 1 and Di is 
fixed at the transversal position or a t 2^,1 ■ We as- 
sume that detectors D s and Di are placed just behind 
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the new slits at the plane-z. These diagonal elements are 
then 



PL s Li,L s Li 



l s lim s mi=zL 



(A) 



u 



(21) 
± for 



where L s ,Li = 0,1 for A, p — x,y or L s ,Li 
X,p = z, and U£" i are coefficients of the effective trans- 
formation Uj given by ([2U)l . One can now obtain the 
non-diagonal elements of the density operator, p, just by 
inverting the above linear equations. Besides, in case of 
a QTR for spatial qudits, the number of necessary slits 
at the plane-z is equal to the dimension D of the qu- 
dits. The positions of these slits can be determined by 
using the eigenvectors of the D 2 — 1 generators of su(D) 
algebra [Uj]. 



A. Diagonal Elements 

The measurement in the logical base, namely p^ z - z \ 
can be determined by coincidence measurements with the 
detectors just behind the double slits [l3[ or at the plane 
of image formation when lenses are placed in the path 
of the double-slits transmitted photons as showed in fig- 
ure QJb) [lil ]. In this last configuration, when the detec- 
tor Dj is at position x z ,q = —100 pm (x Zt i = +100 /an), 
it detects all photons that cross the slit + (— ). By do- 
ing these measurements and normalizing the coincidences 
recorded for the four slits we obtained: = 0.028, 

p+-,+- = 0.468, p-+,-+ = 0.462 and p__ _'_ = 0.042. 



B. Non diagonal Elements 

The second type of coincidence measurements were 
done by positioning the signal detector at x z ,o = —100 
pm or x Zj i — +100 pm. and with the idler detector in 
the plane-z at the transversal positions x y fl = —A/2 
(—0.688mm), x y> \ = A/2 (0.688mm), x x .o = 0mm and 
x Xj i = A (1.376mm) (See figure 0Jc)). When the idler 
detector is at the transverse position x y ,o or x yt i, the 
detector selects the idler photons in the \f(x y o))i state 
or in the|/(a; y> i)) s state. When the idler detector is at 
the transverse position x x fi or x Xt i, it detects the idler 
photons in the \f(x Xj o))i state or in the |/(:Ex,i))s state. 
With these eight measured coincidence numbers, we de- 



termined the non diagonal elements of the density oper- 
ator {/?++,+_, — ,p — ._+}. By repeating 
this detection procedure and reversing the roles of the sig- 
nal and idler detectors (See figureQJd)), we found the non 

diagonal elements {p++. \-,p |_,++,/9_| t --,p--^ }. 

We show below the explicit expressions that determine 
P++-+ 



Re(P++— 



and 



Im(p ++ _ + ) 



where 



Pg+fi+ - P++++ cos 2 



p — i — h sin 9 X 



sin 29 x 



(22) 



-Po+,o+ + P++++ cos2 e v 



p_+_+sin 



sin 26„ 



(23) 



-(Xp,o)\ 



r+(x^ Q )\ 2 + Ir-^i)! 2 



(24) 



The above expressions for Rc(p + ^ |_) and lm(p+-\ |_) 

are obtained by inverting Eq. (|2"Tj) , with L s = 0, Li = +, 
both for A = x, p = z and for A = y, p = z. 

In the third measurement type shown in figure [lje), 
signal and idler detectors are positioned in the detection 
plane-z at the positions x\,k and x^j, with A, p being 
x or y and k,l being or 1. This allows, by means of 
similar expressions to (|22[) and (|23|) . the determination 
of p^ i |_,p }. This set of mea- 

surements correspond to local operations being applied 
to each of the down-converted photons, simultaneously. 



C. The Reconstructed Density Operator 

By performing the quantum tomographic reconstruc- 
tion, as described above, we found the following form for 
the density operator in its matrix representation of our 
experiment 



0.028 0.083 + 0.004i 
0.083 - 0.004i 0.468 
0.081 - 0.005i 0.444 + 0.058i 
-0.129 - 0.062i 0.097 + 0.008i 



0.081 + 0.005i -0.129 + 0.062i ' 

0.444 -0.058i 0.097 - 0.008i 

0.462 0.096 - 0.006i 
0.096 + 0.006i 0.042 



The elements of a density operator must satisfy the H — h, H — , — h and , if it really represents a quantum 

Schwarz inequality, i.e., \pjk\ < y/PjjPkk, where j, k = 
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state. This is not our case for the matrix clement p + ^ , 

since it can be seen that \p + ^ | > ^/p ++++ p The 

reason for that are the experimental fluctuations present 
in the coincidence measurements which can affect the fi- 
nal result as we discussed in the Introduction. This dis- 
crepancy can be reduced by increasing the detection time. 
Even though our reconstructed density matrix presents 
properties which are not fully compatible with the quan- 
tum state description, it is possible to show that it is con- 
sistent with the theory developed in section [Til This is 
done in the next section, where we also show experimen- 
tal evidences of the good quality of our reconstruction. 



IV. DISCUSSION AND CONCLUSION 



The measured density operator shown in (|25|) can be 
approximately written as 



p = 0.87|$) 11 ($|+0.13|$) 22 ($| 
where the states | $ ) are given by 



(26) 



O.O77e l01 |4 
+O.699e i02 | 



O.7O4e i02 |4 
+ 0.099e^ 3 | 



■>, (27) 



and 



0.514| ++) + 0.502e ie | + 
+0.501e l9 | -+>+ 0.483| 



■>, (28) 



with 0i ~ 02 ^ 03 ~ 4.2 and 9 = 0.07. 

However, the possibility to decompose the density op- 
erator, p, in terms of the projectors of a state, |$)i, 
which has a high degree of entanglement and a state, 
| $)2, that is of the form predicted by (0, is not suffi- 
cient for associating them with the states generated by 
each arm of the interferometer in our experiment. We 
still have to give an experimental evidence which corrob- 
orates with (|26p as a reasonable approximation for the 
quantum state of the twin photons, i.e., we need to show 
that the values of A = 0.87 and B = 0.13, obtained 
mathematically, are reasonable for the probabilities of 
generating these states in each arm. 

We measured the values of A and B by blocking one of 
the arms of the interferometer and detecting the trans- 
mitted coincident photons through the signal and idler 
double-slits. A (B) is the ratio between the coincidence 
rate when arm 2 (arm 1) is blocked and the total co- 
incidence rate when both arms are unblocked. From 
this measurement we obtained, A — 0.85 ± 0.03 and 
B = 0.15 ±0.03. 

Another experimental evidence for the high value of 
A can be found in the fourth order interference pat- 
tern recorded (See figure [5]). The interference pattern is 
recorded by using the configuration shown in figure UJe). 
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FIG. 2: Fourth order interference pattern as a function of D s 
position. In (a) , the detector idler was fixed at the transverse 
position x — mm. In (b), it was fixed at the transverse 
position x = 1376 mm. The solid curves were obtained theo- 
retically. 



Fourth order interference pattern as a function of D s po- 
sition was recorded. In (a), the detector idler was fixed at 
the transverse position x = mm. In (b), it was fixed at 
the transverse position x = 1376 mm. The solid curves 
were obtained theoretically [2lJ, with A and B as free 
parameters. Since the state | $)i is almost a maximally 
entangled state we would exp ect to observe conditional 
interference patterns [U |22| when both interferometer 
arms are unblocked. This would not be the case for 
high values of B. The conditionality can be clearly ob- 
served in our interference patterns. The reason for having 
the probability of generating the state from arm 1 much 
higher than the probability of generating the state of arm 
2 is quite simple. The laser beam that cross arm 1 of the 
interferometer is focused at the double slit's plane- za, 
and the spatial correlation of the generated photons is 
such that it is more favorable to the transmission of the 
twin photon through the slits than it is when the pho- 
ton pairs are generated by the pump beam that cross 
arm 2 [12]. These values can be properly manipulated 
by inserting attenuators at the interferometer. 

These experimental observations confirm the good 
quality of the QTR performed on the two photon state 
and allow us to consider the states |$)i and |$)2 as 
good approximations for the states generated by arm 1 
and arm 2 of the interferometer used. Figure [3] shows 
a histogram of the real part of the matrix elements of 
(a) the measured density operator of ([23]) . (b) the den- 
sity operator given by ([2"o]) and, (c) the predicted density 
operator of section [TTJ The agreement between the pre- 
dicted and the measured density operator is good within 
the experimental errors. The largest error for the di- 
agonal elements is only 3.5%. But, for the non-diagonal 
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FIG. 3: Histogram of the real part of the matrix elements 
for (a) the measured density operator, (b) the density oper- 
ator given by ((26} and (c) the predicted density operator of 
Sec. M 



elements the propagated errors reaches 30% for their real 
parts and up to 65% for the imaginary parts. Again we 
remember that these errors can be decreased by increas- 
ing the detection time. 

In conclusion, we have demonstrated that it is possi- 
ble to generate a broad family of mixed states of spa- 
tial qudits by exploring the transverse correlation of the 
down-converted photons. A statistical mixture of spatial 
qubits were used to show the quantum tomographic re- 
construction performed to measure its density operator. 
The process was discussed in details and experimental 
evidences for the good quality of the reconstruction per- 
formed were showed. Even though we had considered 
the state determination only for the case of qubits, it can 
be generalized and performed in a similar way for higher 
dimension systems in a mixed state. The importance 
of this work comes from the possibility of using spatial 
qudits for quantum communications protocols, where it 
requires the ability to characterize them in the presence 
both of noise source and of an undesired user at the quan- 
tum communication channel. 
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